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The success of primordial nucleosynthesis imposes stringent bounds on 
the abundance of gravitational relics. This is particularly true for gravitinos, 
which - for models with gravitationally mediated supersymmetry breaking - 
are expected to have a mass below the TeV scale and thus to decay only af- 
ter nucleosynthesis has concluded. We discuss the nonthermal production of 
gravitinos in models with several chiral fields, to be able to distinguish the 
mechanisms of supersymmetry breaking and inflation. Our explicit calcula- 
tions show that the superpartner of the inflaton can be significantly generated. 
Gravitinos can be produced both at the preheating stage and at the subse- 
quent inflaton/inflatino decay. We verify that this production is well below 
the nucleosynthesis bound, provided the sector responsible for the present su- 
persymmetry breakdown is weakly (for example, only gravitationally) coupled 
to the one which drove inflation. 
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1 Introduction 



There are several issues in astrophysics and cosmology where very high energy scales seem 
to be involved. For instance, consider the observation of cosmic rays with energy above 
the Greisen-Zatsepin-Kuzmin cut-off, or baryogenesis schemes as leptogenesis and GUT 
baryogenesis. The former may be explained with the decay of very massive particles with 
a lifetime comparable or larger than the present age of the Universe; the latter involve very 
massive right handed neutrinos or GUT bosons, whose out-of equilibrium decay generates 
the tiny baryonic asymmetry suggested by primordial nucleosynthesis. In most cases, 
these heavy particles have masses higher than the one of the inflaton, which typically 
ranges in the interval (10 10 — 10 13 ) GeV in the simplest models of "new" and "chaotic" 
inflation. As a consequence, reheating based on a perturbative decay of the inflaton can 
hardly account for their generation. Fortunately, the last decade has witnessed a drastic 
change in the theory of reheating, with the realization that nonperturbative effects can 
play a very significant role. The nonperturbative decay of the inflaton is typically a very 
quick phenomenon, which is expected to be followed by a phase of ordinary reheating 
or thermalization. For this reason, it is commonly denoted as preheating. Preheating 
of bosons is characterized by a very efficient and explosive creation, due to the coherent 
effect of the oscillations of the inflaton field. This allows significant production even when 
single particle decay is kinematically forbidden. It has been successively realized that 
preheating of fermions can also be very efficient despite of the fact that the parametric 
resonance is in this case limited by Pauli blocking. 

One of the main successes of inflation is the generation of a nearly scale invariant 
spectrum of primordial metric fluctuations, as indicated by Cosmic Microwave Back- 
ground and Large Scale Structure measurements. While the nearly scale invariance is 
a natural outcome of the slow roll of the inflaton field, the observed smallness of per- 
turbations typically requires an ad hoc fine tuning of some parameter of the inflationary 
sector. As an example, in the two most standard models of single field chaotic inflation, 
V ((f)) = m;|(/> 2 /2 and V ((f)) = A0 4 /4, one has to fix m\ ~ 10~ 11 Mp for a massive in- 
flaton (Mp here denotes the reduced Planck mass) or A ~ 10~ 13 in the massless inflaton 
case. Supersymmetry may provide a natural framework for the protection of such small 
couplings, and indeed inflationary supersymmetric models have been widely discussed in 
the literature. Since in chaotic and new inflationary schemes the inflaton field acquires 
values comparable or larger than Mp , complete models should include also supergravity 
effects rather than just global supersymmetry. These effects are mostly relevant at very 
early times, but one has to check that they are consistent with the phenomenology of the 
later evolution of the Universe. The most known example of possible difficulties in this 
direction is certainly provided by the gravitino problem. 

The gravitino is necessary present in supergravity, since it is the supersymmetric part- 
ner of the graviton. It has four physical degrees of freedom, two "transverse" and "two" 
longitudinal. Actually, the longitudinal components are present only for broken supersym- 
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metry, and they are provided by the goldstino field, which is a linear combination of the 
partners of the scalars (or gauge bosons) responsible for the supersymmetry breakdown. 
This mechanism is the "fermionic version" of spontaneous breaking of gauge symmetries, 
and indeed it is known as superhiggs mechanism. If supersymmetry has to solve the 
hierarchy problem, the scale of its breaking, and thus the mass of the gravitino (we re- 
strict here to models of gravitational mediated supersymmetry breaking) cannot be too 
far from the electroweak scale, m 3 / 2 ~ 100 GeV — 1 TeV. If not the lightest supersym- 
metric particle, a gravitino with a mass in this range will decay only after nucleosynthesis 
has ended. The decay typically produces an electromagnetic shower, which can alter the 
successful predictions of standard primordial nucleosynthesis. Indeed, not to conflict with 
observations, the very strong limit 

n /2 -^<io- 13 , (1) 

s 

has to be imposed on the gravitino abundance (in standard notation, here denotes 
the gravitino number density, while s the entropy density). 

In the last twenty years, detailed calculations have been performed on the perturbative 
production of gravitinos after the inflaton decay products have thermalized. One finds 
that, for a gravitino with mass in the above range, a rather strong bound on the reheating 
temperature, 

T rh < (lO 9 - 10 10 ) GeV , (2) 

has to be imposed in order to avoid excessive production. More recently, it has been won- 
dered if preheating can also lead to a too large creation of gravitinos. This issue is, how- 
ever, more complicated. Quite generally, preheating effects lead to the non-perturbative 
production of the fermionic partner of the inflaton, the inflatino, and of any other fermion 
which is strongly coupled to the inflaton field. While one easily finds that the transverse 
gravitino component is only weakly (i.e. gravitationally) coupled to the inflaton, and that 
its quanta are hence produced in a negligible amount, the production of the longitudinal 
component is very model dependent. If there is substantial mixing between the inflatino 
and the longitudinal component of the gravitino, the goldstino, preheating may result 
in an overproduction of gravitinos. During inflation and the beginning of reheating, su- 
persymmetry is mainly broken by the inflaton implying a strong correspondence between 
the inflatino and goldstino at this early stage. However, this correspondence does not 
necessarily hold at late times, since supersymmetry may be broken by other fields in the 
true vacuum of the theory. Indeed, this is probably the most typical situation, as it is 
natural to distinguish between inflation and supersymmetry breaking due to the very dif- 
ferent energy scales associated with the two phenomena. In this case, the final gravitino 
abundance can be much smaller than the inflatino one. 

More accurately, the relic abundance of gravitinos will ultimately be related to the 
strength of the coupling between the inflationary sector and the one responsible for the 
present supersymmetry breakdown. The simplest possibility is to consider a model with 
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two fields coupled only gravitationally, which is a simple prototype of more realistic 
schemes of gravitational mediation of supersymmetry breaking in a hidden sector. Even 
in this simplified model, the calculation of nonthermal gravitino production requires sub- 
stantial work. First, one has to develop a new formalism, to be able to clearly define and 
compute the production in systems with several coupled fields. Subsequently, an extended 
numerical investigation has to be carried out to obtain reliable results. This calculation, 
which we have described in details in the works [p], [| , is here summarized in section 0. 

In accordance with the above discussion, the final gravitino production in these schemes 
turns out to be completely negligible also for the longitudinal component. The nonthermal 
production is much more relevant for the fermionic partner of the inflaton, the inflatino, 
which in some cases can be produced with a significant abundance. Thus, gravitinos may 
be overproduced through inflatino decay, if the channel inflatino — > inflaton (or its scalar 
partner) + gravitino is kinematically allowed. Depending on the relative masses of the 
inflaton and inflatino, significant production may be expected instead by the inverse pro- 
cess. This production may be particularly significant if the inflaton sector is coupled only 
gravitationally to matter, since in this case the above decays will have a rate comparable 
to the one generating the thermal bath. Indeed, in such a scenario, the decay channels 
into gravitinos need to be strongly suppressed. In the work || we have discussed the 
possible kinematic suppression of this channel, which translates into a rather strong lower 
bound on the scale of inflation. We have shown that this bound is satisfied by the simple 
single scale (supergravity) inflationary models. This analysis is summarized in section |3|. 

An extensive list of relevant bibliography can be found in |I|, |2|, |3|] . 

2 Gravitino production at preheating 

The system we are considering has the matter content of two superfields <3> and S, with 
superpotential 

W = ^$ 2 + /i 2 (p + S) (3) 

and with minimal Kahler potential G = JC + \n\W\ 2 , /C = $ t $ + S. The field 
(that is, the scalar component of $ ) acts as the inflaton. In the present model, super- 
gravity corrections spoil the flatness of the potential during inflation, so that additional 
contributions must be relevant during inflation. However, we are interested in the dy- 
namics of the system after inflation, when (0) < M p , and supergravity corrections are 
not important. We then assume that the superpotential (0) is valid at this stage, but we 
nonetheless normalize ~ 10 13 GeV, as required by the COBE normalization of the 
CMB fluctuations for the "usual" chaotic inflation. A model where also new inflation is 
obtained was considered in [[|, and it gives the same qualitative results as the one here 
described. 

The superfield S leads to the breaking of supersymmetry in the true vacuum owing 
to its "Polonyi" superpotential. By imposing (3 = [2 — \/3J M p , one can indeed break 
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super symmetry while retaining a vanishing cosmological constant in the true vacuum, 
where the fields s (the scalar component of the superfield S) has expectation value of the 
order of M p . The gravitino mass in the vacuum is of the order n 2 /M p . In order to have 
a gravitino mass of about 100 GeV (that is the expected value for the gravitino mass in 
gravity-mediated supersymmetry breaking models), fi ~ 10 10 GeV is required. 

Right after inflation the field is oscillating about the bottom of its potential with 
frequency proportional to m^. The time dependent expectation value of <fi acts as an 
effective mass for the Polonyi scalar, which has vanishing expectation value at this stage. 
In this initial period the (time dependent) expectation value of is the main source 
of supersymmetry breaking. The amplitude of the oscillations of the field <fi eventually 
decreases, due to the expansion of the Universe, and for times of the order of rri^h the 
Polonyi scalar starts rolling down towards its true minimum and then oscillates about 
it- □ 

The system is thus governed by two time scales. At "early" times, of the order of m^ 1 , 
the only relevant dynamics is the one of the inflaton sector, that is also the main source 
of supersymmetry breaking. At "late" times, much larger than m^ 2 , the system behaves 
as if it was in its true vacuum, and supersymmetry is broken by the Polonyi sector. To 
be more specific, we define the dimensionless parameter fi 2 = /i 2 / (m^ M p ) ~ m 3 / 2 /m^, 
that gives the ratio of the two time scales in the system. If supersymmetry is supposed to 
solve the hierarchy problem, fi 2 should be of the order of 10~ n . Such a small parameter 
implies a very large difference between the two time scales of the problem, which is a 
source of technical difficulties in the numerical computations. As a consequence, we could 
not study the evolution of the system for such a small value of fi 2 . Thus, we kept it as 
a free parameter and we studied how a variation of fi 2 affects the scaling of the relevant 
quantities. 

During the evolution, both the kinetic and potential energies of the two fields con- 
tribute to the supersymmetry breaking. One can define 

(4) 

with rrii = exp (K M- 2 /2) [diW + M~ 2 di K W] . The quantities /< give a "measure" 
of the size of the supersymmetry breaking provided by the F term associated with the 
i— th scalar field. More precisely, we will be interested in the normalized quantities = 
ff I (fi + /I) ) which indicate the relative contribution of the two scalars. It can be easily 
verified |l| that in the initial stages r$ ~ 1 , while r s ~ 1 at late times. The regime of 
equal contribution is around t = fi~ 2 . 

Let us now consider the fermionic content. We denote the fermions of the two chiral 
multiplets by <p (the "inflatino") and s (the "Polonyino" ) . One linear combination of 
them is the goldstino v , while the one orthogonal to v is denoted by T . Initially, v = <p , 

*We neglect here the Polonyi problem associated with the late times oscillations of S . 
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while v = s at late times. In addition, we have the gravitino field, whose longitudinal 
and transverse component are denoted by 9 and by ipj , respectively. The transverse 
component is decoupled from the other fermion fields, and its quanta are produced only 
gravitationally. The longitudinal gravitino component, which in the super-higgs mecha- 
nism is provided by the goldstino v , is however coupled with T . The computation of the 
occupation numbers of the fermions 9 and T is far from trivial and requires a significant 
extension of the existing formalism for nonperturbative production in the one field case. 
The details of the calculation are reported in || . Here we will only present our results. 
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Figure 1: Evolution of the masses of the two fermionic eigenstates. For illustrative por- 
poises we take fi 2 = 10~ 2 . Notice the different normalizations for the two masses. 

As a starting point, one has to diagonalize (at each time) the coupled 9-T system. 
We denote the two fermionic mass eigenstates by ip\ and ip2 ■ In fig. |] we show f\ the 
evolution of their masses for the specific case fi 2 = 10~ 2 . The most relevant information 
which can be constructed by this evolution is very clear: at late times the fields ipi and 
ip2 have, respectively, the mass of the inflatino and of the gravitino field. That is, at 
late times we have the identification ipi = (j) = T , ip2 = 9 (s = v = , being the 
goldstino). This situation is orthogonal to the initial one, when ipi = 9 , ip 2 = s = T 
((p = v = 0). At intermediate times, when supersymmetry is broken by both scalar fields, 
the gravitino is a mixture of ip\ and tJ) 2 ■ To qualitatively appreciate the evolution of the 
gravitino occupation number, we may consider Nq = ri Aq + r 2 N 2 and the orthogonal 
combination = r 2 Nx + r\ N 2 , where are the relative contributions of the two scalars 
to supersymmetry breaking defined above. 

'''From fig. 0, one may be tempted to identify ipi = <fi and tp2 = s . Although this identification 
is rigorous only at the beginning and at the end of the evolution, it can be used for an "intuitive" 
understanding of the system. 
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Figure 2: Evolution of Ng and Nj- for p 2 = 10~ 2 and k = . See the text for details. 

We show in fig. the evolution of Ni , N% , Ng , Nx for modes of comoving momentum 
k = (the scale factor of the Universe is normalized to one at the end of inflation) and 
for p 2 = 10~ 2 . Notice that (by construction) Ng = Ni at early times, while Ng = N2 
at late ones. We also see that ipi is populated on time scales 1 , while ip2 on time 
scales p~ 2 1 . This feature is common for all jl 2 . We remark that the identification 
9 = r\ ip\ + r 2 ip2 should be taken only as a qualitative indication. However, the most 
relevant identification 6 = ip2 at late times is a rigorous one, as should be clear from the 
above discussion. 

We are now ready to present our most important result: the occupation number of 
the two fermionic mass eigenstates at the end of the process. 

As we have said, the realistic case p? = 10 _ 11 is far from our available resources, so 
we kept fi 2 as a free parameter and performed numerical calculations down to fi 2 = 10 -6 . 
The results for fi 2 = 1CT 11 can be clearly extrapolated from the ones we are going to 
present. Moreover, the case fi 2 = can be studied analytically [Q], and it agrees with the 
limit fi — > that one deduces from the numerical results. For p, = , only the inflatino is 
produced. The mass of the Polonyi fermion does instead vanish identically, so no quanta of 
this particle are produced at preheating 0. Notice that p = corresponds to a situation 
with unbroken susy in the vacuum, and it reproduces the models with one single field 
studied so far. We stress that in this case only the inflatino is produced at preheating. 

As we showed in |], 0, the field ipi has occupation number of order one up to k ~ , 
and then rapidly decreases as k~ 4 . The final spectrum is practically independent of 
the value of p~ 2 . These features are easily explained: the eigenstate ipi is associated to 
the inflatino, that is produced by the coherent oscillations of the inflaton; the inflaton 
dynamics occurs on time scales of the order of mT 1 , and is independent of p 2 . 
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Figure 3: Spectrum of gravitinos at late times. 

The spectrum of ip2 - the late time gravitino - is shown in fig. ^ for different values 
of fi 2 , f| and exhibits instead a more significant dependence on fi 2 . One finds that the 
production of these quanta mainly occurs during the first oscillations of the Polonyi field, 
with a typical physical momentum of the order the Polonyi (~ m 3 / 2 ) mass /x 2 /Mp . The 
total number of particles produced in this case is thus an increasing function of fi 2 . One can 
easily realize that the production is expected to be completely negligible in the physically 
relevant case ft 2 = 1CT 11 . 




3 Gravitino production through the inflaton/inflatino 

decay 

In this section we generalize the results discussed above, considering also gravitino pro- 
duction from the decay of inflaton and inflatino. Since we will always consider a super- 
symmetry breaking scale much smaller than the scale of inflation, it is worth noticing that 
the decays of the inflaton and the inflatino are expected to occur nearly simultaneously. 
As we have said, the potential of the inflaton field during inflation is constrained by the 
magnitude of density fluctuations following from CMB temperature measurements. In 
the slow roll regime, one has to impose 

V^ 4 ^ 0.027 e^ 4 M P , (5) 

■'■These spectra are shown at the time t = 10 fL^ 2 m^ 1 . In the // = 10~ 2 — 10~ 4 cases we have 
continued the evolution further, until the spectra stop evolving. In these cases, we have found that the 
spectra very slightly decrease for t > 10 fi" 2 m^ 1 . Thus, the results shown in fig. || give an accurate 
upper bound on the final gravitino abundance. 



7 



where e 2 = M 2 (V /V) 2 /2 1 is one of the two slow-roll parameters (prime denoting 
derivative with respect to 0) and the suffix 60 reminds us that the two quantities have 
to be evaluated when the scales measured by COBE left the horizon, about 60 e-foldings 
before the end of inflation. As stated in the introduction, this relation uniquely fixes the 
scale of the potential in the simplest models of one field inflation and potential V oc n . 
However, models with a much smaller scale than the one reported there and acceptable 
density fluctuations can be constructed. As follows from eq. (H), in models with one single 
field this can be done at the expense of a small e parameter, that is by taking a very flat 
potential during inflation. Such a flat potential may arise more naturally if more scalar 
fields are present, as for example in hybrid inflationary models. 

Due to this freedom, in this section we discuss the production of gravitinos with 
a generic inflationary scale. To be precise, we denote by A 4 the value of the inflaton 
potential at the end of inflation, when the reheating stage begins. Due to the slow motion 
of during inflation, this scale is typically very close to V^o • During reheating, the 
inflaton field oscillates about the minimum O of V, with V (0o) = . As is typical for 
a massive inflaton, we assume that the quadratic term dominates the Taylor expansion 
of V around 0o . We denote by F the amplitude of the inflaton oscillations at the initial 
time t ~ H^ 1 ~ Mp/A 2 . In the present discussion we keep also the value of F as 
a free parameter. As we will now show, the number density of produced gravitinos is 
more closely related to the inflaton mass rather than to the scale A . For a quadratic 
potential and generic values of A and F , one has ~ A 2 /F . Note that the model 
considered in the previous section is characterized by F ~ Mp and ~ 10 13 GeV. 

We discuss here the case in which the inflaton decays only gravitationally. A more 
general analysis where also nongravitational inflaton decays are considered was performed 
in ||. The gravitational decay rate of the inflaton is given by V ~ m^/M 2 , leading to the 

reheating temperature T r h — (T Mp) 1 ^ 2 ~ ^Jm^/Mp . The thermal bound (0) thus simply 
gives 

m <lO 12 GeV. (6) 

The inflaton number density at its decay can be estimated to be — V (0) / at 
t = T<p = T~ 1 . The inflaton "abundance" at the decay time is thus approximatively given 
by 

The abundance of inflatinos produced nonthermally by the inflaton oscillations is also 
easily evaluated, by remembering that inflatinos are produced at preheating up to physical 
momentum , and with a physical number density of the order 



n t ~ 10 2 m\ 



(t = m^)/a(t) 
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The inflatino abundance at t = is thus 



n-^r^io-^^io-^jgj (9) 

Also the abundance of gravitinos produced nonthermally becomes smaller as de- 
creases, so that it is always negligible. Indeed, gravitinos will still be mainly produced 
at the time t ~ with a typical momentum k ~ ^3/2 > independent of the value of 

. However, a lighter inflaton implies a longer lifetime . Since the quantity 113/2/ p^ 4 
decreases in the time interval 777.3/2 < t < , lowering will thus decrease the final 
nonthermal gravitino abundance. 

We see that the inflaton abundance is always higher than the one of inflatinos. Indeed, 
once the bound (|6|) is respected, the inflatino abundance is smaller than the limit 10~ 13 
reported in equation ([]]). Even if all the inflatino quanta produce a gravitino when they 
decay, this would not lead to an overproduction of gravitinos (the inflatino abundance 
can be higher for a massless inflaton, e. g. with potential V ((f)) oc <p n , n > 2 ; see 0, || 
for details). Gravitino production can instead be significant through the — > <j>G decay, 
if kinematically allowed. 

Indeed, assuming that the branching fraction for this decay is 1/N, from eq. (0) we 
get the gravitino abundance 



^ = V* = i^. (io) 

If we now take iV ~ 100 as suggested by the degrees of freedom of the Minimal Super- 
symmetric Standard Model, the bound (0) requires 

m <lO 9 GeV. (11) 

We see that this constraint is much stronger then the one coming from the thermal 
production (0). 

Of course, if it happens that the decay inflaton — > inflatino + gravitino is kinematically 
forbidden (|m^ — mi\ < JTO3/2), then the bound flTT| ) does not hold. | In addition, if the 
scale of supersymmetry breaking is significantly below that of inflation, i.e., /1 <C A, then 
even though the decay <fi — > <p + G is allowed, it will be naturally kinematically suppressed 
with respect to the other decay channels of the inflaton field. This will open an allowed 

§For the specific choice of a present supersymmetry breaking provided by a Polonyi superpotential ([|) 
this decay channel is indeed kinematically forbidden, since one has || 

l m ^ m 4>\ ~ (i ^ TO 3/2 < 771 3/2 . (12) 

Anyhow, such a relation need not be valid in general. 
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window for even in the simplest scenarios. Independent of the details of the decay, 
the rate will always carry a final state momentum suppression factor. The overall decay 
rate can be written as T ~ (l/m ( f ) )\A4\ 2 (p/m^), where \Ai\ is the amplitude for decay and 
the final state momentum suppression factor is 

2p/m = L _ !K^> + H-W) 1/2 ^ (13) 

Thus in models in which m 3 / 2 <C m^, there will be a significant suppression in the produc- 
tion of gravitinos by either inflaton or inflatino decay (note that additional suppression 
may come from the specific form of the amplitude M. as well). 

If we take into account the suppression factor (O) in the bound for the direct gravitino 



production by inflaton decay, and we combine it with the limit coming from the thermal 
production, we find that the inflaton mass has to lay within the interval [|J 

108&v (wgv)^'"^ lor2GeV - (14 > 

Note that the upper bound on comes from avoiding a too quick inflaton decay, with 
a consequent too high reheating temperature. The lower bound is instead due to the fact 
that for too low inflaton mass the kinematical suppression factor ( ]T3"| ) is no longer capable 
of maintaining a safe small branching ratio for the <p — > + G channel. 

The simplest models of single scale inflation with purely gravitational decays and a 
scale set by (|5[), have typically F ~ Mp and A ~ few x 10~ 4 M P . In this case the 
gravitino mass is slightly higher than 10 10 GeV, well within the allowed range (|I4"D . We 
thus conclude that these models do not suffer from a gravitino problem. 
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